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In the radiation analysis of the NERVA Nuclear Rocket System, two-dimensional 
discrete ordinates calculations are sufficient to provide detail in the pressure 
vessel and reactor assembly. Other parts of the system, however, require three- 
dimensional Monte Carlo analyses. To use these two methods in a single analysis, 
a means of coupling was developed whereby the results of a discrete ordinates cal- 
culation can be used to produce source data for a Monte Carlo calculation. Several 
techniques for producing source detail were investigated. Results of calculations 
on the NERVA system are compared and limitations and advantages of the coupling 
techniques discussed. 


In the analysis of the NERVA** nuclear rocket 
engine, a complex radiation transport problem must 
be solved. Since any shield design must simultan- 
eously minimize weight and radiation level, both 
penetration and geometry effects are important. 
ANALYSIS METHODS 

Two rigorous methods have found wide appli- 
cation in complex radiation transport problems : 
discrete ordinates (Sjj) and Monte Carlo. A wide 
range of computer codes is available for both 
methods. 

The discrete ordinates (Spj) method can be used 
in one -dimensional or two-dimensional problems. It 
employs an iterative balance of radiation in a 
large number of mesh regions. Angular and energy 
dependence is calculated in a limited number of pre- 
selected directions and energy groups. Collision 
angle- to-angle and energy-to-energy transfers are 
determined by expansions of the scattering proba- 
bility in terms of Legendre polynomials. Where 
these approximations can be employed, the method 
allows a relatively rapid evaluation of radiation 
transport throughout a system. If the system con- 
tains voids or regions of widely varying collision 
densities, two-dimensional analyses may result in 
anomalies called ray effects^ in which the 
allowed directions do not properly intercept regions 
of interest. No three-dimensional representation is 
possible. 


*Public Release Approval: PRA/SA - SNPO-C, 
dated 2k November 1970. 

**The Nuclear Engine for Rocket Vehicle Application 
(NERVA) program is administered by the Space 
Nuclear Systems Office, a joint office of the 
USAEC and NASA. Aerojet Nuclear Systems Company 
is prime contractor for the engine system and 
Westinghouse Electric Corporation is principal 
subcontractor responsible for the nuclear sub- 
system. 


The Monte Carlo method may allow pointwlse 
representation of radiation distribution in space, 
energy, and angle with exact kinematics to the 
degree to which cross section data are available. 

Any amount of detail must be paid for in computer 
time. The Monte Carlo method calculates and 
tabulates the fates of individual particles (or 
individual ’'particles" which may represent any 
number of real particles or energy quanta). If 
relatively few collisions occur or averages over 
large regions, angles, or energies are allowed, 
the method can be quite efficient. 

APPLICATION TO NERVA 

A NERVA propulsion module consists of a pro- 
pellant tank and engine which includes pressure 
vessel and reactor assembly, shielding, nozzle, and 
the necessary pumps, pipes, valves, instrumentation, 
and controls. A typical layout is shown in Figure 1. 

The pressure vessel and reactor assembly 
( FVARA ) is the primary source of radiation. For 
radiation analysis puiposes, it consists of a core, 
beryllium reflector, forward internal shield, and 
pressure vessel. The PVARA is essentially a two- 
dimensional object and can be accurately modeled 
in cylindrical geometry. The PVARA materials are 
sufficiently dense so that many mean free paths 
would be seen by particles before leaking from the 
outer surfaces of the pressure vessel. Pure analog 
Monte Carlo analyses of the PVARA would be very 
costly in computer time, especially if statisti- 
cally valid radiation leakage data are needed in 
narrow bands of space, angle, or energy. 

The portion of the engine between the PVARA 
and the propellant tank includes several thousand 
pounds of metallic components in complicated three- 
dimensional geometries with many void regions. In 
the nearly empty condition, the propellant tank and 
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Figure 1 

NUCLEAR PROPULSION MODULE CONFIGURATION 


liquid hydrogen impose few mean free paths to most 
radiation. 

A comparison of the capacities of the analyti- 
cal methods and the characteristics of the system 
indicates that Monte Carlo is necessary for most 
regions but Sjj may be sufficient for analysis of 
the pressure vessel and reactor assembly. 

METHOD COUPLING 

If S N were to be used for analysis of the FVARA 
and Monte Carlo for most other regions of the sys- 
tem, a means must be devised for producing source 
information for Monte Carlo from Sjj results. This 
information must be in sufficient detail for the 
needs of any desired analysis. 

In each spatial mesh in each energy group, an 
angular flux can be calculated by along several 
predetermined vectors. Each such vector can be 
visualized as a line from the center of an igloo to 
the center of a block of the igloo. Figure 2 indi- 
cates how such an igloo would describe the direc- 
tions in which the angular flux could be described 


in a single radial band on the top surface of the 
reactor cylinder. The Sjj fluxes are tabulated in 
energy, space, and angle. One could define a 
"particle” for each such vector in terms of the 
particles /second value of the tabulated flux. These 
vectors, if continued by Monte Carlo through a void 
or low density region, would not necessarily inter- 
sect regions of interest. This is comparable to 
the "ray effects" obtained if the S N model contains 
voids. To more nearly represent the physical prob- 
lem, the particle directions must be distributed 
uniformly over that part of the igloo that can be 
seen from the regions of interest. 




Figure 2 

S N ANGULAR DESCRIPTION 

HISTOGRAM 

The easiest way that particles can be produced 
to emerge between the vectors predicted by Sjj is to 
sample a uniform distribution of vectors and assign 
to each vector that emerges through any block its 
share of the particles/sec value of the Sjq- vector 
in that block. This is equivalent to a histogram 
in angle with constant particle current within the 
solid angle intercepting the block. Tracing the 
particles evaluated in this mode results in a 
spatially oscillating flux at a distant surface. 

In particular, radiation coming from the end of the 
PVARA will generally be underpredicted near the 
polar axis and radiation from the side will be 
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overpredicted directly above the outer radius of 
the reactor. 

SMOOTHING 

In order to improve on the histogram results, 
a system of least-square curve fits in polar and 
azimuthal angle was devised. This results in a 
smoothly varying assignment of particles/sec to 
vectors at all angles and produces a smoothly 
varying flux at distant surfaces. The values at 
the polar axis and above the reactor radius, in 
general, are improved by the smoothing procedure 
but are limited by the extrapolation of the curve 
fits from the nearest vector. 

TAILORED QUADRATURES 

The quadrature set, the set of angles defining 
the vectors, can be tailored to a particular 
analysis to remove the need for extrapolating curve 
fits over angles of several degrees. This can 
improve the accuracy by providing that several of 
the "blocks" of a more refined igloo can be seen 
from the detector surface, 

COMPARISON OF METHODS 

A computer program was written to accept sur- 
face angular flux data from the Sjj code, D0 t( 2 ) 

and produce the necessary fitting and sampling. 

(U) 

This code, DASH W ’ ' can also calculate the trans- 
port of the sampled vectors through a void or pure 
absorber. 

As an example of the results from histogram, 
smoothed, and tailored quadrature modes, a series of 
calculations was performed using an S^ quadrature 
set and an quadrature set*. The results from 

the S£ calculation were treated in the histogram 
mode and the smoothing mode. Since the re- 

sults were tailored to use very small steps in 
angle in the direction of interest, only the histo- 
gram mode was needed. The Sg quadrature used a 
smallest polar angle of 22 degrees while the small- 
est angle was 1.4 degrees. 

-X- 

The quadrature order given here is that con- 
ventionally used and is not entirely consistent. 

The Sg set contains 24 vectors distributed 
symmetrically and the set contains 12 

symmetrically distributed vectors directed down- 
ward and 9k vectors directed upward with many 
near the polar axis. 


Figure 3* shows the dose rate in rad(c)/hr 
from gamma radiation leaking from the top of the 
FVARA to a plane 37 meters above it. A void is 
assumed in the 37 meters to accentuate the dif- 
ferences between the results of the different 
methods. It is seen that the S£ histogram method 
predicts values well below those of the S6 smoothed 
and which agree very well for this case. The 

smallest S£ angle intercepts this plane at about 
1600 cm. from the axis. The histogram mode, then, 
includes the assumption that the particle current 
over the plotted range is the same as that ex- 
pected at a much larger radius. The smallest Sqp4 
angle, on the other hand, intercepts this plane at 
about 90 cm * Since the source position has a range 
given by the 140 cm diameter of the reactor, any 
such underpredictions will be masked by compensat- 
ing overpredictions. Since the S5 smoothed results 
are produced by a curve fit which provides an ex- 
trapolation from the vector at 1600 cm, the excel- 
lent results are perhaps fortuitous. Poorer 
smoothed results would be expected when the angular 
flux values of the different vectors varied 
more than in this case. 
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Figure 3 

GAMMA DOSE RATE FROM PVARA END AT PLANE 
IN VOID 37 METERS ABOVE PVARA 

*While these data are obtained by collecting par- 
ticles in radial interval bins, they are shown 
here by smooth curves in order to make comparisons 
easier. 
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For a more difficult case, the radiation is 
considered which leaks from the side of the reactor 
and moves upward at a small angle with the side. 

The same three calculational modes were employed 
and transported through a void to a plane 3 *^ 
meters above the reactor. This is approximately 
the location of the bottom of the propellant tank. 
These results are shown in Figure 4 . 

The detector plane in this calculation, at 
tank bottom, is close enough to the source loca- 
tions to be intercepted by the vectors at the 
first two polar angles of the Sg quadrature set. 

As a result, the oscillatory behavior of the 
histogram mode is exhibited and two peaks are seen 
which correspond to the two angles. 



GAMMA DOSE RATE FROM PVARA SIDE AT PLANE 
IN VOID 3. A METERS ABOVE PVARA 

Considering only those particles which can 
strike the bottom of the propellant tank, smoothing 
does not significantly improve the S£ data. This 
is seen to be due to the difficulty of using data 
from larger angles to extrapolate to angles smaller 
than 22 ° from the side of the reactor. The side 
leakage at small angles can be considered further 
by extending the particles to the plane at 37 meters 
as was done for the top leakage. These results are 
shown in Figure 5. 

As in the case shown in Figure 4 , these re- 
sults are for angles smaller than the vector 
nearest the polar axis and: give an indication of 
the accuracy of the use of Sg quadratures for a 
difficult application. 
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RADIAL POSITION, cm 

Figure 5 

GAMMA DOSE RATE FROM PVARA SIDE AT PLANE 
IN VOID 37 METERS ABOVE PVARA 

COUPLING INTERNAL TO THE Sfl PROBLEM 

In considering the detail of small angle re- 
sults from an Sj _24 calculation, it appears that 
the intended detail of the problem may be extreme. 
Some scatter was observed in flux versus angle plots 
from the S^ij. calculation. With regard to the 
approximations inherent to an calculation, such 
detail probably is not completely reliable. Hie 
loss of strict kinematic angle-energy relationship 
in using expansions in Legendre polynomials causes 
S$ calculations to predict some results in violation 
of energy-angle restrictions. As an example of 
this, the side leakage case was considered further 
in order to discover the origin of small angle 
radiation. The outer 2 centimeters of the PVARA is 
the aluminum pressure vessel wall. Because of the 
steep angle at which radiation from inside the wall 
must traverse the aluminum in order to be directed 
toward the plane at 37 meters, the low energy 
gammas are strongly attenuated. 

To study the effect of the wall, Sjj angular 
fluxes from the DOT code using an S5 quadrature 
set were considered in the radial interval just 
inside the wall. The DASH code was then used to 
accomplish two calculations. First, DASH was used 
to sample particles in the direction of the 37 
meter plane, attenuate them through the wall, and 
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calculate their transport to the plane. This gives 
the dose rates at the plane from radiation which 
does not scatter in the wall. Second, DASH sampled 
particles in all outward directions into the wall 
and wrote these data on magnetic tape suitable for 
input to the COHORT Monte Carlo code. The wall was 
modeled in COHORT and the tape used as a source 
to calculate the transport of scattered gammas to 
the 37 meter plane. The results of these two 
calculations were summed* and compared with the S 5 
smoothed and S 124 calculations using angular fluxes 
at the outer surface of the wall. The energy 
spectra in a band at 400 - 500 cm on the plane are 
shown in Figure 6. The Monte Carlo result compares 
well with the S^24 result over most of the energy 
range. Apparently the attenuation of the wall at 
very small angles cannot be properly accounted for 
in the S£ calculation, resulting in a higher pre- 
dicted dose rate. Probably the decoupling of angle- 
energy relationships in the procedure hinders 
the assignment of scattered radiation to the lowest 
energy group, accounting for the steep drop in dose 
rate near zero energy. The peak at 6 to 7 Mev is 
due to sources in the reflector of the reactor. 
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S 124 FROM SURFACE 


MONTE CARLO FROM 
S 6 INSIDE WALL 


E, Mev 


Figure 6 


ENERGY DISTRIBUTION OF GAMMA RADIATION 
FROM PVARA SIDE IN 400-500 CM RADIAL BAND 
ON PLANE IN VOID 37 METERS ABOVE PVARA CENTER 


*The effect of gamma sources in the wall was con- 

sidered also, but did not contribute significantly 
to the total. ^ 


In a similar application^ sin S£ DOT calcula- 
tion was used to produce angular flux information 
between the reactor core and reflector. The par- 
ticles simulated at this interface were traced 
through the reflector and wall regions by COHORT 
to a surface above the reactor > and a data tape 
of particle information produced. A pure Monte 
Carlo calculation of sources in the reflector and 
wall was also run using COHORT to produce a second 
data tape of particle information above the reactor. 
These data tapes were then used with a COHORT model 
of a propellant tank to calculate radiation levels 
at the top of the tank with different quantities of 
liquid hydrogen in the tank. 

A similar set of calculations was performed 
using S£ smoothed data at the outer surface of the 
reactor to produce a tape of particle information 
above the reactor. This tape was then used by 
COHORT to calculate radiation levels at the top of 
the tank with several different quantities of 
liquid hydrogen propellant in the tank. The results 
are compared for a particular hydrogen level in 
Figure 7. In this case, the tank has a capacity of 
about 180,000 lbs and is about 1/3 full. The top 
of the tank is about 30 meters from the top of the 
reactor. The results are divided into that radi- 
ation which does not scatter in the tank structure 
or contents (direct) and that which does scatter 
(scattered). These dose rates are consistent with 
earlier results in indicating the slight over- 
prediction inherent in the data. Since leakage 
from the side of the reactor is the primary contri- 
butor to dose rate at the tank top, it is not sur- 
prising that the relationship between surface 
source and core-reflector interface source results 

should be similar to the relationship between the 
S 5 and S-j g 4 data seen earlier. 

OTHER EXAMPLES OF THE METHOD OF CODE COUPLING 

The use of the method is not limited, as we 
have seen, to using surface angular fluxes. Nor is 
it limi ted to tracing radiation through a void* 

Source data can be produced at any grid within the 
Sjj problem, axially or radially. The data can be 
produced in the particle source format used by 
many Monte Carlo codes, (it has been used with TMC, 
COHORT and MORSE. ) It can also produce tab ula r 
source data for FASTER or additional DOT problems. 
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Figure 7 

CAM4A DOSE RATE AT TANK TOP FROM PVARA SIDS 
SOURCES — BASED ON DOT CALCULATION 


ADVANTAGES OF THE METHOD OF CODE COUPLING 

The method allows considerable flexibility in 
choice of techniques to solve many complex radiation 
transport problems. It allows a user to take ad- 
vantage of many of the special characteristics of 
Sjj and Monte Carlo codes. Where the detail of 
Monte Carlo is needed, it can be employed. Where 
is sufficient, it 'can be used with significant 
saving in computer time. Different degrees of 
detail are available in by using tailored quad- 
rature sets. 

Any number of samples can be produced at a 
problem interface. As an example of the advantage 
of this, consider a problem of the type in which 
one determines the effect of scattering in the wall 
on the radiation level 37 meters away. A large 
number of particles scattering in the wall is neces- 
sary since few of the scattered particles reach 
detector regions offering small solid angles to the 
scattering region. If the particles must originate 
in a pure Monte Carlo sense, where the real parti- 
cles originate, few ever arrive at the scattering 


region of interest. Thus, a great many original 
particles would have to be produced in order to 
scatter a few in the wall. 

But many scatters are necessary in the wall. 
Probably this problem would be prohibitively costly 
in computer time for all but the most liberally 
financed establishments. Using the Sjj - Monte 
Carlo coupling technique, however, the problem can 
be quite economical. A low quadrature can be 
run on an IBM 360/65 in an hour or less for a 10 - 
20 energy group, 50 x 30 spatial mesh problem. A 
high quadrature problem, such as an S-j_ 2 ^, may take 
two to three times as long. The DASH problem to 
produce Monte Carlo data or calculate transport 
through a void or pure absorber takes 5 to 10 
minutes. DASH can produce any required number of 
source particles. A Monte Carlo problem in the 
wall would take a few minutes longer. The whole 
process takes little more time than the original 
S N . Furthermore, once the % problem has been run 
and the angular flux data saved on magnetic tape, 
it can be used any number of times. This allows 
additional analyses tailored to specific needs at 
any later time without repeating the bulk of the 
calculations. later calculations can be performed 
in any degree of Monte Carlo employment, from pure 
Sjj to pure Monte Carlo, with commensurate computer 
costs. 

The coupling method is of general applicability 
and can be used to produce useful source data for 
any Monte Carlo analysis. It is most efficient for 
calculating radiation environment from shielded or 
distributed sources. Typical applications would 
be in calculating fluxes or dose rates in a void, 
air, or complicated three-dimensional structures 
from reactors or shielded isotopes. 
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